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To understand the seemingly absent temperature dependence in the conductance of two- 
dimensional topological insulator edge states, we perform a numerical study which identifies the 
quantitative influence of the combined effect of dephasing and elastic scattering in charge puddles 
close to the edges. We show that this mechanism may be responsible for the experimental signatures 
in HgTe/CdTe quantum wells if the puddles in the samples are large and weakly coupled to the 
sample edges. We propose experiments on artificial puddles which allow to verify this hypothesis 
and to extract the real dephasing time scale using our predictions. In addition, we present a new 
method to include the effect of dephasing in wave-packet-time-evolution algorithms. 


The discovery of two-dimensional topological insula¬ 
tors, 2d-TIs, stirred up hope for potential applications 
using their ability to host reflectionless one-dimensional 
transport. In practice, however, the transport along 
the edges of the state-of-the-art realizations of 2d-TIs 
exhibits a length-dependent resistance which increases 
above the expected quantized value in samples signifi¬ 
cantly larger than 1 /rm. Moreover, most notably, the 
edge state resistance seems to exhibit no observable 
temperature dependence in a wide temperature range 
(T « 30 mK — 30 K) [1]. This seems to be hard to rec¬ 
oncile with the notion that the backscattering on 2d-TI 
edges is due to inelastic processes. Still, this behavior 
has been been frequently observed in various recent ex¬ 
periments on 2d-TIs, both based on HgTe/CdTe [2-4] 
and InAs/GaSb samples [1, 5]. Understanding this ap¬ 
parently generic feature is important for the field, as it 
might show up ways to improve the edge transport in 
existing material systems. 

As coherent elastic backscattering is symmetry forbid¬ 
den at 2d-TI edges [6], there have been theory studies 
considering alternative backscattering mechanisms. Be¬ 
sides the possibility for elastic backscattering by mag¬ 
netic impurities which explicitly break time-reversal sym¬ 
metry [7-9], they mainly focused on inelastic backscat¬ 
tering by Coulomb interaction or phonons on the edge 
[10-18] as well as Coulomb backscattering in quantum 
dots which are expected to appear naturally along the 
edges of the currently used material systems (HgTe/CdTe 
and InAs/GaSb quantum wells) due to trapped charges 
at the gate insulator interface [19, 20]. 

In addition, there have been a few proposals in which 
the inelastic interaction is mainly held responsible for 
breaking the electron phase coherence while the backscat¬ 
tering is then attributed to elastic scattering. For a de¬ 
phasing process which does not explicitly flip the spin, it 
was found that one does not expect backscattering along 
a clean edge of a 2d-TI ribbon, as backscattering requires 
a full spin flip in this setup [21]. This changes if a quan¬ 
tum dot—in which the extended states are naturally spin 
mixed due to the spin-orbit coupling—is present along 
or close to the edge. In this setup, even a decoherence 


mechanism which conserves the average spin may lead to 
backscattering, as has been qualitatively shown in Ref. 22 
where dephasing was mimicked by adding lattice sites 
with an imaginary self energy. 

In this publication, we also address dephasing-induced 
backscattering. However, given the fact that the in¬ 
teresting T-independent edge state resistance has been 
consistently observed in various systems and material 
classes, we do not attempt to assign it to a specific 
source for decoherence along with a corresponding spe¬ 
cific T-dependence of the dephasing/decoherence time. 
Instead, assuming that edge channels are coupled to 
nearby puddles due to spatial charge inhomogeneities 
[3, 19, 20, 22, 23] we present a more general study which 
is aiming at a quantitative description of backscattering 
through the interplay of dephasing and spin mixing in 
such quantum dots. This includes to consider puddle 
dwell times and couplings to the puddles, which are cal¬ 
culated for realistic model Hamiltonians for the case of 
the HgTe/GdTe material system. Thereby, we extract a 
range of dephasing times that would be compatible with 
the measured temperature-independent resistance, and 
which is a function of the puddle dwell time. This de¬ 
pendence of the conductance on the puddle properties 
could be checked in future experiments with artificially 
created puddles. 



FIG. 1. (color online) Snapshots of the calculated wavefunc- 
tion (the color encodes the spin) in the device geometry. A 
wave packet is approaching a puddle (top panel) in which it 
is fully spin-mixed by spin-orbit scattering (bottom). A video 
of the dynamics including dephasing is available online [38] . 
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FIG. 2. (color online) Time-dependent probability density in a pnddle of length L close to the edge which is traversed by a 
wave packet starting on the edge close to the puddle at t — Ops, see Fig. 1. (a) Behavior for short times showing the steep 
rise in density from the wave packet entering the puddle and a small fraction which is leaving the puddle after traversing it 
ballistically. Most of the density decays linearly (black dashed fit lines), for longer times shown in (b). The inset shows the 
extracted linear decay coefficients run (see main text) from the fit as a function of L together with a power-law fit. (c) Log-log 
plot of the derivative of the density (the current) demonstrating also the power-law decay of the density for long times. 


For our calculations, we use the Bernevig-Hughes- 
Zhang (BHZ) Hamiltonian [24] discretized on a square 
grid. Intrinsic spin-orbit coupling of the Dresselhaus type 
is included, which mixes the two spin blocks [25]. We 
study the time evolution of a wave packet which is ini¬ 
tially localized on the edge of a 2d-TI and approaches 
a nearby puddle defined by a local electrostatic poten¬ 
tial close to the edge, cf. Fig 1, employing a method 
which has been shown to be well suited to describe the 
edge-state dynamics in HgTe/CdTe 2d-TIs [26, 27]. The 
detailed calculation setup is described in detail in the 
Appendix. 

From the time evolution of the wave packet, we cal¬ 
culate the time-dependent probability density p{t) to 
find the charge carrier in the puddle, which includes 
edge-puddle-coupling and lifetime effects. Impurity- 
configuration-averaged results of such calculations for 
various puddle sizes are shown in Fig. 2. As can be seen 
from the short-time dynamics. Fig. 2(a), the wave packet 
enters the cavity on a time scale of a picosecond. A small 
fraction of the density, 2 -12 % depending on the puddle 
size, exits the puddle after a ballistic traversal. However, 
most of the density stays in the puddle and decays only 
slowly with an almost constant outflow, visible as a linear 
density decay for intermediate times, cf. Fig. 2(b). The 
shown density is scaled to the total density of the wave 
packet and the fact that values close to one are reached 
reflects the forbidden backreflection on the clean edge: 
Most of the density has to couple into the puddle in this 
geometry. 

As soon as the absolute density in the puddle falls be¬ 
low Ri 0.5, the time dependence of the density changes 
into a power law. This is best seen in Fig. 2(c), which 
shows the negative time derivative of the density in a log- 


log plot. Here, one clearly recognizes the linear density 
decay as a plateau with superimposed oscillations, which 
can possibly be attributed to ballistic orbits in the pud¬ 
dle. This turns over into at~~ -power-law decay for long 
times as indicated by the dashed lines [28]. Such power- 
law behavior is known [29] and similar results have al¬ 
ready been obtained in wave-packet simulations on other 
material systems [30]. 


With only a few additional assumptions, the knowledge 
of the dwell time of the electrons in the quantum dot will 
allow us to estimate the effects of dephasing on its trans¬ 
port properties. The first one is that the decoherence dis¬ 
turbs the electron phase evolution but does not explicitly 
lead to spin flips, which is quite realistic for setups with¬ 
out magnetic fields and a very low density of magnetic 
impurities. Due to the strict spin-momentum locking of 
the TI edge states, this spin conservation implies that 
the dephasing process will not lead to backscattering as 
long as the electrons are propagating along the edges. 
In a quantum dot however, the spin is fully mixed al¬ 
ready after a very short time, due to the combined effect 
of impurity scattering and intrinsic spin-orbit coupling 
(see Fig. 1). Thus, the efficiency of backscattering due 
to dephasing will depend on the relation of the dephas¬ 
ing timescale and the electron lifetime in the puddle. To 
make this quantitative, we assume that decoherence is 
statistically occuring with independent events. We first 
additionally assume that each event leads to full phase 
loss (an assumption that we will later refine). Thus, the 
electron density being in the puddle at the event will 
leave it in a random direction in the subsequent dynam- 

























3 


ics. Then, the average reflection can be estimated by 
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where is the mean time between dephasing events and 
^max is the time at which the density pit) in the puddle 
attains its maximum pmax- Equation (1) makes use of 
the piecewise monotonous structure of p{t) and can be 
understood as follows: Suppose that the first dephasing 
event occurs at a time t 2 > tmax with the density already 
decaying monotonously. As each event is assumed to 
lead to full phase randomization, all subsequent events 
will not matter as the propagation in the puddle is al¬ 
ready fully random. Hence, the knowledge of the first 
event after reaching pmax is enough to determine the to¬ 
tal backscattering probability for the time window of the 
decay. It can be calculated as an expectation value of 
pit) with an exponential distribution for the dephasing 
event, describing the mean waiting time in a Poisson pro¬ 
cess. This main contribution enters Eq. (1) as the second 
term. The first term in Eq. (1) describes the backscat¬ 
tering during the (monotonous) rise of pit). Here, the 
argument can be reversed as the latest event will deter¬ 
mine the amount of backscattered density in the time 
window of the rise. The last term in Eq. (1) takes care 
of the double counting that occurs for events both on the 
rising and the falling edge of pit). 

For the shape of the density curves obtained from the 
numerical calculations, one can simplify the above ex¬ 
pression. Since the rising edge is very short compared 
to the decay, it suffices to consider only the second term 
of Eq. (1). In addition, we showed above that the sub¬ 
sequent decay is approximately linear. The power law 
decay can be omitted as long times are exponentially 
suppressed in the integral. We then find for the aver¬ 
age transmission, using pit) ss pmax [1 - (t - tmax)/Tiin], 
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J 
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Here rnn is the time after which the puddle would be 
empty assuming a steady linear decay. The values for 
T[in which were extracted from the linear fits to pit) are 
plotted in a log-log plot in the inset of Fig. 2(b). One 
finds that Tiin approximately scales like 

Using these formulae, we can estimate the average 
transmission probability for a single puddle of size L for 
a given decoherence time as shown in Fig. 3. As antic¬ 
ipated, we find good agreement between the results from 
the full model, Eq. (1), and the approximation with a 
linear fit, Eq. (2). Both show a characteristic sigmoidal 
behavior with a saturation at T = 1 — /Omax /2 for small 
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FIG. 3. (color online) Average transmission of a single quan¬ 
tum dot at the edge as a function of the dephasing time for 
different dot sizes. The full lines are evaluated from the full 
time-dependent density, using Eq. (1), while the dotted lines 
show the result from Eq. (2), assuming a linear decay. 


T 0 . For the smaller puddles, the extracted value for pmax 
is underestimated for the linear fit, see Fig. 2(a). 

More generally, the value of Pmax is not only given by 
the puddle size, but also by the coupling of the puddle 
to the edge and it would therefore be strongly depen¬ 
dent on the distance of the puddle from the edge. This 
coupling will also influence the dephasing-time scale be¬ 
low which one observes the nearly constant backscatter¬ 
ing probability: With worse coupling, there will be less 
backscattering per puddle but the saturation will already 
be reached earlier. The extracted cutoff times therefore 
refer to the “perfect coupling limit” meaning that they 
may be interpreted as a lower bound for the cutoff time. 

So far, the dephasing only entered the model a posteri¬ 
ori but was not included in the calculation of the electron 
dynamics. We found a way to include it in the wave- 
packet dynamics calculation using an algorithm which is 
inspired by the concept of einselection [31]. It treats de¬ 
phasing in a rather general fashion and allows not only to 
vary the strength of single dephasing events but also to 
study the influence of dephasing on the wave-packet dy¬ 
namics itself. For details of the implementation, we refer 
to the Appendix. The underlying interaction responsible 
for the dephasing is left unspecified and the dephasing 
time remains as a parameter. Also, the implementation 
may not be applicable for all kinds of environments but 
we think that it captures the important effects and ex¬ 
pect similar results for other implementations. 

Our results for the total transmission of dynamical de¬ 
phasing calculations are shown as symbols in Fig. 4(d). 
Here, the dephasing was chosen to be less effective such 
that one needs on average 2.75 events to achieve full de¬ 
phasing (the time scale of the dephasing time axis corre¬ 
sponds to full dephasing). The dotted lines show a fit to 
the data using Eq. (2) which expectedly agrees well for 
small dephasing times. In this limit, the time between 
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FIG. 4. (color online) (a)-(c) Time-dependent density in the 
puddle obtained including different effective dephasing times 
in the time evolution. The colors correspond to different pud¬ 
dle geometries and couplings as shown in the inset in (d). (d) 
Averaged transmission through the puddle as a function of 
the chosen dephasing time for the different puddle geometries. 
The dashed lines are fits to the numerical data (symbols) us¬ 
ing Eq. (2). 

the events is small compared to the dynamics of the sys¬ 
tem, thus many small events lead to the same net result 
as one strong event. For longer dephasing times, how¬ 
ever, in our refined calculation using partial dephasing, 
it comes into play that the dynamics of the system will 
be strongly influenced by the dephasing as this effectively 
opens another exit channel for the puddle thus decreas¬ 
ing the electron lifetime. This can be nicely seen in the 
plots of the time-dependent density for different dephas¬ 
ing times in Fig. 4(a~c). In the long dephasing time limit, 
many frequent weak events lead to more backscattering 
than one rare strong event, as the second exit channel 
(backscattering) will be (partially) opened already at an 
earlier time. 

Fig. 4 also shows the influence of different puddle ge¬ 
ometries on the transmission and the system dynamics. 
The results match the expectation that puddles further 
away from the edge will lead to less backscattering in 
the strong dephasing limit but will also saturate below a 
higher dephasing time cutoff. 

With this, we can now relate to the experimental re¬ 
sults on 2d-TI samples. So far, the measured resistance 
does not show any observable temperature dependence in 
the experimentally accessible T-range (30mK-30K) [1]. 
This would be in line with our findings if the dephasing 
times in the experimental samples were below the cutoff 
time even at the lowest temperature. Then, an increase 
of the temperature (and a connected reduction of the 
dephasing time) would not lead to additional backscat¬ 
tering. To make a specific example, if the backscatter¬ 


ing in the sample was mainly caused by 500 nm puddles 
which are well coupled to the edges, this would imply 
that the dephasing time should be shorter than « 10 ps. 
For a rough comparison, from low-temperature magne¬ 
toconductance measurements on 2d HgTe quantum well 
samples, dephasing times of ~ 70 ps were extracted [32]. 

Thus, given that the coupling to the puddles in the real 
samples is likely to be smaller, it might indeed be that 
the lack of observed temperature dependence is mainly 
due to the inherently short dephasing times in these ma¬ 
terials. We only consider puddles of a depth of 40meV 
as we are limited to the range of validity of the BHZ- 
Hamiltonian. Deeper puddles, similar to larger pud¬ 
dles, would also shift the curves toward longer dephas¬ 
ing times and decrease the temperature above which the 
conductance is expected to be temperature independent. 
The puddle depth and the connected change of electron 
density may also influence t^, however, at least in the 
2 d-limit, recent experiments show controversial results 
whether an increasing density leads to an increased [32] 
or a decreased [33] dephasing time. Note that, in line 
with experiments, our model predicts reproducible con¬ 
ductance fluctuations which are similar to UCFs as the 
coupling to the puddles is expected to be depending on 
gate voltage and magnetic field. However, the conduc¬ 
tance is not T-dependent as long as the dephasing time 
is below the cutoff. 

To make a definite statement, one would have to ex¬ 
perimentally check this hypothesis which could be done 
with experiments on artificially created puddles. As there 
are experimental samples which show the quantization, 
it seems to be possible to create puddle-free samples in 
which one could artificially introduce puddles of varying 
size using a local gate. This should allow for experi¬ 
mentally reproducing the calculated sigmoidal curve and 
extract the temperature-dependent dephasing time. 

To summarize, we quantitatively studied the backscat¬ 
tering of 2d-TI edge states due to the interplay of dephas¬ 
ing and dynamical scattering in charge puddles. Our re¬ 
sults suggest that the seemingly absent temperature de¬ 
pendence could be due to the saturation of backscattering 
at dephasing times which are short compared to the pud¬ 
dle lifetimes. One should be able to verify this hypothesis 
quantitatively with experiments on small (« 100 nm) ar¬ 
tificial puddles, which, according to our study, should 
show a detectable temperature dependence at the ex¬ 
perimentally available temperatures and from which one 
could extract the actual dephasing time. 

In passing, we developed a scheme to implement 
dephasing into wave-packet-time-evolution algorithms 
which is general and could be used in a wide range of sce¬ 
narios. It is described in detail in the Appendix. For the 
charge puddles in 2d-TIs, it would also allow for the ex¬ 
plicit inclusion of other sources of backscattering, like ex¬ 
ternal magnetic fields, which we recently showed to have 
a strong influence on the backscattering due to puddles 
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Appendix 


Model Setup 


In the main text, we show numerical calculations of an 
edge-state wave packet which is scattered by an electro¬ 
static puddle and extract the time-dependent probability 
density p{t) of the wave packet being localized in the pud¬ 
dle. We model the electronic structure of a HgTe/CdTe 
quantum well using the Bernevig-Hughes-Zhang Hamil¬ 
tonian [24], 


H = 


( 


h{k) 

0 zi 


\ -Z\ 0 


0 -z\ \ 

z\ 0 

h*{-k) ^ 


(3) 
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with spin-subblock Hamiltonians 
J 


Mk) = 


V(x) + M(x) - {B + D)k^ 
Ak- 


Ak^ 

V{x) — M{x) -\- {B — D)\^ 


(4) 


where k± = kx i iky and = k'^ + ky. Intrinsic 
spin-orbit coupling of the Dresselhaus-type is included 
by the parameter A in Eq. (3), which mixes the two spin 
blocks [25]. We use an electrostatic potential V{x) to 
model circular and stadium shaped puddles and confine 
the states by the potential M[x) to get a quantum spin 
Hall edge state at the system boundary. For the calcula¬ 
tions presented in the manuscript, the potential strength 
of the puddle was set to 40 meV leading to hole-like 
states. 

As initial state, we create a Gaussian edge-state wave 
packet which is assembled in reciprocal space along 

the boundary. Throughout the manuscript, we use a 
Gaussian with a width in position space of cr = 90 nm. 
Also, we include only states propagating towards the 
electrostatic puddle, leading to a strongly spin-polarized 
wave packet. We calculate the propagation of the wave- 
function using a propagator based on an expansion 
of the time-evolution operator in Ghebyshev Polynomi¬ 
als [35]. During the propagation of the wave packet tpit), 
we integrate over the puddle region resulting in 

the time-dependent probability density p{t). In order to 
avoid effects of mesoscopic fluctuations, we average over 
20 different configurations, which differ in a random im¬ 
purity potential with an amplitude of 5 meV and a wall 
distortion of 20 nm. 

Propagation with dephasing 

In the following section, we will show how dephas¬ 
ing is included in the numerical calculations presented 
in Fig. 4. The dephasing algorithm described here is 
inspired by the concept of einselection (“environment- 
induced superselection”) pioneered by W. Zurek [31], 
which is a mechanism proposed to understand the influ¬ 
ence of dephasing on quantum systems and, more gener¬ 
ally, to explain the quantum-to-classical transition. Ac¬ 
cording to einselection, the interaction of an open system 
with its environment leads to decoherence [36], which 
causes a decay of the quantum states of the system into 
an incoherent mixture of so-called pointer states. This 
strongly suppresses quantum interference effects between 
different pointer states on the time scale of the dephas¬ 
ing time T 0 . The character of the set of pointer states, 
the pointer basis, heavily depends on the environment 
and the coupling to it. For example, for very weak cou¬ 
pling to the environment, the pointer basis coincides with 


the set of energy eigenstates of the system. However, in 
the case of an intermediate system-environment coupling 
based on a local interaction, e. g., in the case of electron- 
phonon or electron-electron coupling, it is expected to be 
a set of states that is localized in phase space, i. e., in 
position and momentum. 

As in the puddle lifetime calculations, we again employ 
a numerical time evolution based on a single pure state, 
i. e., we do not use a representation in terms of density 
matrices, which would drastically increase the computa¬ 
tional effort. Still, we incorporate the dephasing-induced 
interference suppression by occasionally (partially) ran¬ 
domizing the phases of the components of the state vector 
in a representation that tries to faithfully mimic a decom¬ 
position in terms of the pointer basis. This randomiza¬ 
tion is done at fixed event times which are sampled 
from an exponential distribution with time constant Te, 
i. e., we assume the events to be fully uncorrelated (Pois¬ 
son process). In between these events, the propagation 
is done fully coherently using the propagator based on a 
polynomial expansion mentioned above. The decompo¬ 
sition and the subsequent randomization is done in the 
following way: At the time of an event we extract a 
set of pseudo eigenstates, 

nt„+At 

OC / '0(t) (5) 

Jtn-At 

at the energies Em = {0 meV, ±1.5 meV,... ± 7.5 meV} 
using a short-time propagation of the wave packet ipii) 
around the time tn- These states fulfill the requirement 
that they are both local in energy as well as in position 
space (as they are extracted from a propagation over a 
finite time interval). The degree of localization in posi¬ 
tion space can by tuned by the propagation time At [37], 
which we choose At ~ 1 ps. We use these states to 
spectrally decompose the current state before the event 
ip{tn), as sketched in Fig. 5. For this, we first orthogo- 
nalize them with a Gram-Schmidt process, leading to the 
set (pm- Then, we calculate the weights am^ which can be 
used to express the current state as the decomposition 

= '^amilm+Alp. (6) 

m 

Since the set of 11 states is not sufficient to describe the 
wave packet tp{tn) exactly, we also allow for account a 
small residual part Aip. In this decomposed representa¬ 
tion, the dephasing can be easily added by modifying the 
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FIG. 5. The figure illustrates the workings of the dephasing 
algorithm. At a dephasing event, the current wavefunction— 
here showing a wave packet spread in an electrostatic puddle 
(the color encodes the spin)—is spectrally decomposed in a set 
of pseudo eigenstates and recomposed with random phases, 
leading to a new dephased wavefunction. A video using this 
type of dephasing is available online [.38]. 


set of amplitudes to 

am = am exp(i7r rand), (7) 

with the random numbers rand S [—Q,Q]- Subsequently, 
we create the new wave packet 

+ e) = X! (8) 

m 


which is used in the following propagation. With this 
dephasing algorithm, we are able to calculate the time 
evolution in arbitrary mesoscopic systems. It roughly 
conserves the energy, the spatial extent, as well as the 
total spin of the wave packet, which is what one would 
expect from the interaction with a spin-unpolarized (non¬ 
magnetic) environment. The strength of a single dephas¬ 
ing event can be tuned by the degree of randomization 
of the phase in each event, i. e., by the parameter Q. 
This also makes the regime of weak but frequent dephas¬ 
ing accessible (compared to the regime of rare and very 
strong dephasing which allows a simpler model treatment 


discussed in the first part of the manuscript). However, 
to compare with the experimentally accessible dephasing 
time Tcj), which is the time for full phase coherence loss, 
one has to find the relation between the time constant of 
the events and the time for full dephasing . This can 
be done in the following way: The average correlation of 
a state in the pointer basis with itself after n events 
is given by, 



■ dqr^ 


E?=i ij 


1 

(^ 


Q 


dqi ■ ■ ■ dqn cos 




'-Q 


dq cos {irq) 


/ sin (ttQ) \ 

I ) 


(9) 


With the chance for n events occurring after time t in a 
Poisson process. 


P{n) 


r 1 __L 

--e 

n\ 


( 10 ) 


one can evaluate the average decay of autocorrelation 
after time t. 


(^P(0)|^/>P (t)) coherent 


n—0 

oo 


E l / t sin {ttQ) 

n I 


n—0 

= exp 

= exp 


ni \Te ttQ 
t sin (ttQ) t 

Te 


_Te TtQ 


t 

Te. 


1 - 


sin (ttQ) 
ttQ 


= exp 


t 


( 11 ) 


with 


T0 = 


ttQ 


ttQ — sin (ttQ) 


Te, 


( 12 ) 


yielding the desired relation between decoherence time 
T 0 and the time constant of the events Tq. To model the 
frequent but weak dephasing regime, we use Q = 1, thus, 
T0 = 7r/(7r - 2)re « 2.75 Te. 

In Fig. 4, we apply the dephasing scheme on three dif¬ 
ferent electrostatic puddle shapes to show how the cou¬ 
pling between the puddle and the edge states affects the 
total spin randomization of the puddle in presence of de¬ 
phasing. A video of a sample propagation can be also 
found online [38]. 





















